Non-equilibrium dynamics in a 3d spin-glass 
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Non-equilibrium dynamics in a Ag(ll%Mn) spin-glass has been studied by low frequency ac- 
susceptibility and magnetic relaxation experiments. The results unequivocally show that spin struc- 
tures that memorize the cooling process are imprinted in the system. These imprinted structures 
disclose themselves through dramatic changes of the dynamics on re-heating the spin-glass through 
the temperatures where intermittent stops or changes of the cooling rate have been imposed. We 
can qualitatively interpret our results in terms of the droplet spin-glass model developed by Fisher 
and Huse [Phys. Rev. B 38 (1988) 373; 386]. 
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I. INTRODUCTION 

A successful real space theory for the dynamics of 
3d spin-glass systems is the droplet model developed a 
decade ago by Fisher and HuseB Particular experimental 
characteristics such as the aging phenomenonB and the. 
non-existence of a phase transition in a magnetic fieldEI 
are inherent properties of the droplet spin-glass model. 
However, the validity of this model has been questioned: 
e.g. a compact domain growth picture is not consistent 
with the multiple jncmory effect recently found in 3d 
spin-glass systems.cllj Fig. 1 shows a low frequency ac- 
susceptibility experiment on a Cu(Mn) spin-glass where 
this suggestive effect is illustrated. In this paper, ex- 
tended measurements of some associated non-equilibrium 
features of the ac-susceptibility and the dc-magnetic re- 
laxation on the 3d spin-glass material Ag(llat% Mn) are 
presented. The results can qualitatively be incorporated 
in the original droplet modelEJ without apparent contra- 
dictions and we do not here try to interpret our results on. 
aging in spin-glasses in terms of non real space models 



II. THE DROPLET MODEL 

The results of this study are interpreted in the spirit 
of the droplet model using relevant parts as summarized 
in the following points; 

i) The spin-glass ground state is unique but two-fold 
degenerate by its spin reversal symmetric state. 

ii) Chaos with temperature - a small temperature shift 
changes the ground state spin configuration com- 
pletely on long enough length scales.Q 

iii) The length scale, I, up to which no essential change 
in spin configuration is observed after a tempera- 
ture change AT* introduces the concept of an over- 
lap length /(AT*). 

iv) In the non-equilibrium case, the development to- 
wards the ground state is governed by the growth 



of domains belonging to either of the two degener- 
ate ground states. The typical domain size after a 
time, t w , at a constant temperature T, is 



R(x 



T\n(t w /T ) 
A(T) 



(1) 



-13 



where tq is a microscopic time of order 10~ iu s. 
A(T) sets the free energy scale and ip is a bar- 
rier exponent. This length scale R, that we often 
will refer to as a typical domain size, should not be 
taken literally, but rather be considered as a typical 
measure of the smallest separation between domain 
walls remaining after the wait time t w . 

In order to discuss the relaxation towards equilibrium, 
a 3d Ising spin-glass quenched to a temperature T\ < T g 
is considered. The state at t = is hence characterized 
by spins of random direction. To further simplify the fol- 
lowing argumentation it is also assumed that the spins 
are situated on a regular lattice. In the droplet model 
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FIG. 1. The out-of-phase component of the ac susceptibil- 
ity of a Cu(13.5 at% Mn) (T g w 68 K) spin-glass measured on 
continuous heating after the sample has been subjected to two 
intermittent stops of order 10 4 s at 40 and 50 K during the 
cooling procedure. The two pronounced dips in the curve oc- 
cur at the temperatures of the stops at constant temperature 
during cooling. / = 1.7 Hz. 
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there exists only two degenerate ground states, & and 
its spin reversal counterpart 9, it is thus possible to 
map all spins to either of the two desired ground states, 
immediately after a quench. Citing Fisher and HuscEJ 
"This results in an interpenetrating network of regions 
of the two states". In fact, the domain walls of one 
lattice spacing, which separates the two states <P and 
<F at t = after the quench should have a fractal sur- 
face giving rise to fractal domains of all sizes. The sub- 
sequent equilibration process at t > is governed by 
droplet excitations yielding a domain growth to a typi- 
cal length scale which depends on temperature and wait 
time, t w i, according to Eq. [jj After this time, fractal 
structures typically smaller than R(T\,t w i) have become 
equilibrated, but, viewing the system on a larger length 
scale, it still looks intact. The implication of this argu- 
ment is that large enough fractal structures which were 
obtained at t = after the quench are virtually unaf- 
fected after the wait time t w ±. If the system thereafter is 
quenched to a considerably lower temperature T2, a new 
fractal domain structure can be mapped onto the equi- 
librium configuration at T 2 . The short length scale limit 
is now set by the overlap length, l(Ti — T2), and frac- 
tal structures will be washed out with increasing wait 
time, t W 2, starting from this length scale and ending at 
the size R(T2,t W 2)- Structures on longer length scales 
are again intact. If the system is further quenched and 
then instantly re-heated to T2, the domain structure at 
T2 of course remains unaffected, but also, when instan- 
taneously heated back to T\ , the re-structuring on short 
length scales, R(T 2 ,t W 2) <C R(Ti,t w x), that occurred at 
T2 is rapidly washed out and the system is effectively left 
with only the original large length scale domain struc- 
ture which was imprinted at the first quench to T\. If 
the fractal droplet model is relevant, these basic prop- 
erties of the logarithmically slow (eternal) equilibration 
process should be reflected in the dynamic response of 
a real 3d spin-glass subjected to a corresponding ther- 
mal history. In a real experimental situation, the system 
can only be cooled and heated at a controlled but finite 
rate. The domain growth is in such processes governed 
by the cooling/heating rate and limited by an interplay 
between the domain growth (Eq. [j]), the chaotic nature 
of the spin glass phase and the overlap between states 
on short enough length scales. Accounting for this com- 
plication, we have designed a series of dynamic magneti- 
zation and susceptibility experiments that should mirror 
the discussed development of a domain structure in an 
equilibrating 3d spin-glass. 

Droplet excitations occur on time limited length scales 
and are in fact responsible for the equilibration process 
of the quenched spin-glass, also giving rise to the do- 
main growth function Eq. [I]. If a weak magnetic field 
is applied on the spin-glass at time t b s = 0, the time 
dependent response, m(t obs ), is due to a continuous mag- 
netization process governed by polarization of droplets of 
size, L(t obs ): 



L oc 



Tln(t obs /T ) 
A(T) 



(2) 



Since L grows with the same logarithmic rate as R, the 
relevant droplet excitations and the actual domain size 
become comparably large at time scales t w w t f, s . For 
t bs *C t w the relevant excitations occur mainly within 
equilibrated regions while for t obs 3> t w these excita- 
tions occur on length scales of the order of the grow- 
ing domain size and yield a non-equilibrium response. A 
crossover occurs in the intermediate region, t w ~ t obs 
which is characterized by a maximum in the relaxation 
rate S{t) = 1/H dm(t)/d\n(t). This aging behavior is 
also reflected in low frequency ac-susceptibility measure- 
ments. Different from dc-relaxation measurements is that 
the observation time is kept constant by the probing fre- 
quency, u/27r, according to t obs — 1/lo and hence, in an 
experiment at constant temperature, sets the associated 
probed length scale, L(T,l/u), fixed. In this study we 
are primarily interested in the processes that affects the 
response at the observation time of the experiment, the 
consequences of these are best exposed in plots of the 
relaxation rate S vs. t obs for different wait times, the 
out-of-phase component of the ac-susceptibility vs. time 
at constant temperature or vs. temperature at different 
cooling/heating rates. Therefore, the dc-relaxation mea- 
surements are primarily presented in the form of the re- 
laxation rate, S, and the ac-susceptibility is represented 
by its imaginary part—x", these quantities can be related 
through the relation,!! 
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x" ~ 2 S ^ obs ) ' tobs = 1 / UJ 



III. EXPERIMENTAL 
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The sample was prepared by melting pure Ag and Mn 
together at T = 1000 °C in an evacuated atmosphere. Af- 
ter annealing the sample at 850 °C for 72h it was water 
quenched to room temperature. 

The experiments have been performed m a non- 
commercial low-field SQUID magnetometer The dc 
magnetic field is generated by a small superconducting 
solenoid always working in persistent mode during mea- 
surements. The ac-field was generated by a copper coil 
directly wound on the sample which is shaped as a 5 
mm long cylinder, 2.5 mm in diameter. The magnetic 
response of the sample, subtracted by the ac-field from 
a compensating coil, was recorded with a set of pick-up 
loops positioned to form a third order gradiometer. At 
the position of the sample, the resulting rms value of 
the ac-field was 0.1 Oe and the background field was less 
than 1 mOe. All ac-susceptibility measurements were 
performed at a frequency of 1.7 Hz. A sapphire rod 
was used to provide a good thermal contact between the 
heater, the thermometer and the sample. 
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IV. RESULTS AND DISCUSSION 



A. Aging characteristics 



Figs. 2 and 3 introduce the overall behavior of the 
dynamic susceptibility of our Ag(llat%Mn) sample and 
expose two classical manifestations of the aging phe- 
nomenon in spin-glasses. In the inset to Fig. 2 both com- 
ponents of the ac-susceptibility, x( w ) = x'( w ) +x"( UJ ) at 
the frequency iojlix = 1.7 Hz, are plotted vs. tempera- 
ture at a cooling rate of 0.25 K/min. The maximum in 
X'{T), defining the freezing temperature, is accompanied 
by the onset of an out-of phase component of the suscep- 
tibility. In the main frame of Fig. 2 the influence of aging 
on x" is illustrated by a plot of \" vs. time elapsed at 23 
K, where the continuous cooling was interrupted. Fig. 3 
exemplifies the results from an ordinary zero field cooled 
(ZFC) aging experiment at T m = 27 K, in (a) the mag- 
netization and in (b) the corresponding relaxation rate is 
plotted vs. the observation time. The sample is probed 
in a dc-field of 1 Oe after the wait times t w = 100 s (open 
symbols) and t w = 3 000 s (filled symbols) at T m . The 
signatures of aging in spin-glasses: an inflection point in 
m(t) and a corresponding maximum in S(t) vs. \ogt b s 
are reproduced in the figure. 

Some words about time and non-equilibrium dynam- 
ics: experimentally we are confined to time scales set 
by the cooling/heating rates and the time it takes to 
thermally equilibrate the sample to a constant temper- 
ature. Although spin-glass dynamics occurs on all time 
scales ranging from 10~ 13 s to infinity, the shortest ob- 
servation times where consequences of non-equilibrium 
dynamics can be observed experimentally are confined 
to t bs > 10~ 3 t ae ff where t ae ff is an effective age of the 
system set by the cooling/heating rate or the time spent 



at constant temperature. I.e. when t a b s <C taeff no ag- 
ing is observed and the dynamic response of the system 
appears stationary. In practice this applies to most ac- 
susceptibility measurements at frequencies larger than 20 
Hz (t b s = 1/u) « 10~ 2 s). On the other hand, in a dc 
magnetic relaxation experiment spanning some decades 
in time, the observation time of the experiment con- 
tinuously increases and log t b s unavoidably approaches 
logt ae ff at long enough observation times and effects 
of the aging process becomes discernible. The different 
times necessary to discuss the aging behavior are related 
according to: t ae ff « t w + t Q b s = t a , where t a is the total 
time spent at constant temperature i.e. the age of the 
system. 

We will continue by discussing results from ac- 
susceptibility experiments using different cooling rates 
and employing intermittent stops at different tempera- 
tures during the cooling of the sample, and interpret this 
cooling behavior as well as the subsequent heating curves 
(always recorded at one and the same heating rate) in 
terms of an assumed domain structure imprinted during 
cooling. This discussion will be complemented by results 
from dc-magnetic relaxation experiments subsequent to 
cooling and heating procedures that closely mimic those 
of the ac-experiments. 
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FIG. 2. Relaxation of x"(t) with time at constant tem- 
perature after the sample has been cooled from a tempera- 
ture above T g to T — 23 K. Inset, both components of the 
ac-susceptibility vs. T measured at a cooling rate of 0.25 
K/min. / = 1.7 Hz, h ac = 0.1 Oe, Ag(ll at% Mn). 
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FIG. 3. Zero field cooled magnetization vs. observation 
time at 23 K. a) m(t) vs. \ogt b) relaxation rate S(t). H dc = 1 
Oe, Ag(ll at% Mn). 
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B. ac-susceptibillity: cooling rate dependence and 
memory 

A simple memory effect in the low frequency ac- 
susceptibility is shown in Fig. 4 after the following se- 
quence has been performed. The sample is cooled to 23 
K while measuring every 0.25 K (filled circles). At 23 
K the sample is subjected to a wait time t w — 7 200 s 
during which x" decays (cf. Fig. 3). After this wait 
time, the sample is further cooled to 20 K and then im- 
mediately re-heated, the data on heating are indicated 
by open circles in the figure. Both the cooling and the 
heating is made with the same rate 0.25 K/min. The 
results from identical cooling and a heating procedures, 
but without the wait time at 23 K, are indicated by the 
full lines in Fig. 4 for reference. The effect of the re- 
laxation during the wait time is clearly visible in the 
cooling curve as a dip in x"(T), followed by a rather 
rapid increase and approach to the reference curve at 
temperatures below 23 K on the continued cooling and, 
as expected from the droplet model, it also reappears as 
a corresponding dip in x"(T) centered around T = 23 K 
in the heating curve. The merging with the cooling refer- 
ence curve on continued cooling.-can be interpreted as an 
effect of chaos with temperature!!! and a region of overlap- 
ping states, AT*. The dip on re-heating is a consequence 
of the imprinted spin configuration at T — 23 K on large 
length scales. I.e. when returning to this temperature 
all domain growth at lower temperatures have occurred 
on length scales <C i?(23 K,7 200 s) and these are here 
washed out on shorter time scales leaving a system with 
an effectively equivalent domain configuration to the one 
originally obtained during the wait time at T — 23 K. 

An immediate consequence of the chaos and overlap 
concepts is that the ac-susceptibility depends on the 
cooling rate. Larger domains have time to grow within 
the overlap region AT* when cooling the spin-glass at a 
slower rate. Thus, the observed magnitude of x"(T) will 
be closer to the equilibrium level in a slow cooling process 
than in a fast. This prediction is confirmed in the inset 
to Fig. 4, where x"(T) at two different cooling rates, 
0.25 and 0.005 K/min. are shown. The overlap con- 
cept is further supported by the fact that when suddenly 
at 25 K the cooling rate is increased from 0.005 to 0.25 
K/min., the cooling curve within a distance of AT* w 1 
K merges with the fast cooling curve. Assuming that the 
domain size, R(T,t c ), at all temperatures in the cooling 
process is limited by a characteristic time t c , which is set 
by the cooling rate and an interplay between chaos and 
the overlap length, obviously this time t c must increase 
with decreasing cooling rate. In the x"(T) experiments, 
droplet excitations on the length scale L(T, 1/w) govern 
the response, these excitations includes a large amount on 
the size of the domains if 1 /to is of order t c , but such ex- 
citation rapidly becomes rare with an increasing t c . This 
implies as mentioned above that the measured suscepti- 
bility decreases with decreasing cooling rate and that 
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FIG. 4. x"(T) vs. T measured on cooling (solid circles) and 
heating (open circles) at a rate of 0.25 K/min. The sample 
was intermittently kept at T — 23 K for t w = 7 200 s during 
cooling. Solid lines show the corresponding curves without 
the interrupted cooling at 23 K. The inset shows x"(T) using 
two different cooling rates: 0.25 K/min. (open circles) and 
0.005 K/min. (solid circles). The slow cooling rate is changed 
to 0.25 K/min. at 25 K. / = 1.7 Hz, h ac = 0.1 Oe, Ag(ll 
at% Mn). 

a change of the cooling rate must result in a crossover 
from one cooling rate characteristic curve to the other. 
It is of course also possible to cross from the fast to the 
slow cooling curve in a similar way. 

Regarding the two reference curves, measured at 0.25 
K/min., in the main frame of Fig. 4, it—may be ar- 
gued that they should be almost identical^ since a do- 
main structure on continuously decreasing long length 
scales R(T, t c ) is imprinted in the system on cooling the 
sample. This structure should then also govern the ac- 
susceptibility on re-heating the sample at the same rate. 
However, as is seen in the figure, the two curves differ 
notably in magnitude. Two main factors can explain the 
deviation. Firstly, the effective heating time scale, th, 
is different from the cooling time t c , although the rates 
are the same, 0.25 K/min. This difference can be un- 
derstood considering the existence of an overlap length 
and the fact that the domain growth, or, equivalently, 
the domain wall movements are faster at higher tempera- 
tures. At any subsequent temperature, T, on a cooling or 
heating process, overlapping ground state configurations 
have already been created at nearby temperatures within 
T ± AT* . Hence, the domains have grown to a larger size 
when reaching T in cooling and a smaller size on heat- 
ing, which in terms of time implies: t c > th- If the map- 
ping was perfect, so that the domain structure on long 
length scales was fully intact on re-heating, the heating 
curve would still coincide with the cooling curve as long 
as th <t c , since we would probe an identical undisturbed 
domain structure determined by R(T,t c ) (> R(T,th)) at 
all temperatures. However, secondly, there is a reinitial- 
isation of the gained domain structure on length scales 
that do not overlap with the created domains at a lower 
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FIG. 5. x"(T) measured on heating at a rate of 0.25 
K/min. The different curves are recorded after cooling 
the sample at different cooling rates: 0.005 K/min. (solid 
squares), 0.25 K/min. (open circles) and 4 K/min. (solid 
circles). / = 1.7 Hz, h ac = 0.1 Oe, Ag(ll at% Mn). 

temperature in a continuous cooling process. This im- 
plies that the domain structure on re-heating the sample 
has become partially reconstructed also on length scales 
of order R(T, t c ), i.e. the system appears less equilibrated 
and the susceptibility attains a comparably larger magni- 
tude. The behavior on changing from cooling to heating 
the sample requires a comment. In Fig. 4 this change 
occurs at T — 20 K resulting in a heating reference curve 
which initially is smaller in magnitude than the cooling 
reference curve. This is explained by the fact that the 
sample spends the longest time in the temperature inter- 
val within the overlap region (20 K +AT*) ± AT* first 
during cooling and then in the subsequent heating. 

To elaborate further on the cooling/heating rate de- 
pendence, we show in Fig. 5 three curves obtained at one 
and the same heating rate, but recorded after cooling the 
sample at 4 K/min., 0.25 K/min. and 0.005 K/min. The 
differences between the curves are significant and in ac- 
cord with the discussion above that a domain structure 
on long length scales R(T,t c ) is imprinted but also reini- 
tialized during continuous cooling. We thus expect that 
the curve recorded at a cooling rate of 4K/min. should 
always yield th > t c and the system should appear least 
equilibrated and have the largest magnitude of %"(T). 
The curve obtained after cooling at 0.005 K/min. has 
an imprinted spin structure characterized by t c 3> th 
and x"(T) should have a substantially smaller magni- 
tude. The curve measured after cooling at 0.25 K/min. 
is expected to be found in-between the other two, just as 
is seen in the figure. 

For visual clarity in the following figures, all presented 
ac-susceptibility data are subtracted by either the cool- 
ing or the heating reference curve. These two reference 
curves are obtained using a cooling/heating rate of 0.25 
K/min. (see Fig. 4). Since the ac-susceptibility always 
decreases towards the equilibrium value when the system 
ages at constant temperature, data points which are 
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FIG. 6. x"{T)-x"e.f{T) vs. T, measured at a cooling (solid 
symbols) or heating rate (open symbols) of 0.25 K/min. The 
sample was kept 30 000 s at 24 and at 21 K during cooling. 
The open squares indicate an experiment where the sample 
only was kept 30 000 s at 24 K during cooling. / = 1.7 Hz, 
h ac = 0.1 Oe, Ag(ll at% Mn). 

positive in such plots indicate a state further away from 
equilibrium than the corresponding points on the refer- 
ence curve. For negative values the data points are closer 
to equilibrium. Furthermore, in all these experiments, 
the de-field is zero and the cooling is always initiated at 
T = 40 K where the system is paramagnetic. 

As has been experimentally observed in a recent study 
on memory effects in spin-glassesB and was reproduced in 
Fig. 1, the single memory experiment presented in Fig. 
4 can be extended to include two (or more) wait times 
at well separated temperatures during the cooling of the 
sample. A procedure, which in the subsequent heating 
process results in dips in the ac-susceptibility positioned 
at each of these temperatures. For our current sample, 
this intriguing behavior is illustrated in Fig. 6 where the 
sample has been subjected to two wait times of 30 000 
s duration, first at T = 24 K and then at T = 21 K. 
Except at these particular temperatures the sample is 
cooled and heated with the usual rate of 0.25 K/min. As 
can be seen from the open circles in Fig. 6, the equili- 
bration at each temperature gives rise to corresponding 
dips when re-heating the sample. 

Evidently from the results discussed above, it is pos- 
sible to imprint a spin configuration at a temperature 
where the spin-glass has been allowed to age and recover 
this state when returning to the same temperature, and 
even to imprint and recover two or more equilibrated 
spin configurations at well separated temperatures. This 
ability excludes, as pointed out recentlyEl that new com- 
pact domains grow at each aging temperature. In other 
words, assuming that \T\ — T 2 \ > AT*, it is impossible to 
imagine that compact T\ domains coexist with compact 
T2 domains. However, by assuming that the equilibration 
process does not correspond to the growth of compact do- 
mains but rather, as discussed above, a removal of fractal 
domain wall structures up to the length scale R, the 
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FIG. 7. x"(T) - Xref(T) vs. T, measured on cooling at 
different rates (solid symbols) and heating (open symbols) at 
a rate of 0.25 K/min. Solid circles indicate a cooling process 
at 0.005 K/min. down to 25 K, where the cooling rate is 
changed to 0.25 K/min. Open circles show the corresponding 
curve on heating at 0.25 K/min. Solid squares below 25 K, 
show data after having kept the sample at 25 K for 7 200 s us- 
ing a cooling rate of 0.25 K/min. otherwise and open squares 
the corresponding heating curve. / = 1.7 Hz, h ac — 0.1 Oe, 
Ag(ll at% Mn). 

results may still accord with the droplet model. To be 
more specific, the minimum length scale of the distance 
between domain walls, after the first wait time^at T± is 
typically R(Ti,t w \). Chaos with temperatureQ implies 
that this structure is random at T% on length scales larger 
than the overlap length l(T\ —T2). After the wait time, 
t W 2, at T2, the domain wall movements have removed 
fractal structures on length scales up to R(T2,t W 2)- The 
spin configuration of the system can now be mapped 
on the equilibrium configurations at T2 and at T\. The 
map at T2 yields an equilibrium system on length scales 
smaller than R(T2,t W 2) and a random system on larger 
length scales. The map at T% on the other hand shows 
random structures on length scales between l(Ti— T2) and 
R(T2,t W 2) and on length scales larger than R(Ti,t w i). 
In the real experiment, the configurations are less well 
defined and how well they are preserved in the ther- 
mal procedure depends on the length of the wait times, 
the cooling/heating rates and the separation between the 
temperatures T\ and T2. On re- heating the sample, the 
relative depth of the recovered dips in x"(T) in Fig. 6 
mirrors how these structures have been preserved. For 
comparison, the ac-susceptibility on heating the sample 
from 20 K after waiting t w = 30 000 s only at T = 24 K 
during the cooling, is included as squares in Fig. 6. Since 
this single memory dip is deeper than the corresponding 
double memory dip it can be concluded that during the 
wait time at T = 21 K the first imprinted spin configura- 
tion at the preceding aging at T — 24 K has been partly 
over written, i.e. the system has been partly reinitial- 
ized. It is also worth to notice that the relative depth of 
the dip on re-heating the sample is temperature depen- 



dent; with increasing temperature, in otherwise identical 
single memory experiments, the dip becomes more and 
more shallow compared to the depth obtained during the 
original wait time at T m on cooling (see also Fig. 7). 

A partial reinitialisation, is also evidenced from the 
positive values of the relative ac-susceptibility plotted in 
Fig. 6 a few Kelvin above and below the dip at 24 K, 
signaling an apparently less equilibrated system as com- 
pared to the reference. Before discussing the origin of 
this behavior more thoroughly, it is advantageous to re- 
turn to the very slow cooling experiment. In Fig. 7 
the ac-susceptibility during cooling with 0.005 K/min., 
subtracted by the corresponding data obtained with the 
cooling rate 0.25 K/min., is represented by filled circles. 
When passing the temperature T — 25 K the cooling 
rate is changed to 0.25 K/min. resulting in a quick in- 
crease of the susceptibility towards the reference level. 
In fact it even crosses the reference curve during the pro- 
ceeding cooling to T = 20 K. The ac-susceptibility of 
the subsequent heating, at rate 0.25 K/min., arc shown 
by the open circles. As can be expected when consid- 
ering the overlap length and a decrease of the reinitial- 
isation just when the cooling rate is increased, the sys- 
tem looks most equilibrated in a temperature interval 
T = (25 K+AT*) ± AT*. At higher temperatures, the 
spin configuration appears less equilibrated than the orig- 
inal slow cooling curve, but still closer to equilibrium 
than the reference heating curve (which is measured at 
the same heating rate, 0.25 K/min.). This behavior mir- 
rors both the continuously imprinted equilibrated spin 
configurations obtained during slow cooling and the si- 
multaneous partial reinitialisation of these configurations 
which occurs at lower temperatures. It must hence also 
be concluded that the equilibrated length scales R(T, t c ) 
obtained at each temperature in the slow cooling process 
are partly preserved during the whole experiment. For 
comparison we have in Fig. 7 included a single memory 
curve recorded at a cooling/heating rate of 0.25 K/min. 
and waiting 7 200 s at 25 K during cooling. Filled squares 
show the cooling curve and open squares the heating 
curve. It may be noticed that this dip is centered around 
T = 25 K as expected. 

Fig. 6 and 7 show that a long wait time at constant 
temperature also affects the ac-response at surrounding 
temperatures, so that the spin-glass appears less equili- 
brated than in the reference procedure at temperatures 
separated more than AT* from the temperature where 
the aging took place. As discussed above, in the con- 
tinuous cooling case the typical equilibrated domain size 
R(T, t c ) is governed by the effective cooling time t C) which 
also sets an effective age and thus the magnitude of 
x"(T,t c ). When the cooling is stopped at T m , the do- 
mains can grow unrestricted, and after the wait time 
they have reached a size R(T m ,t w ). When cooling is 
recaptured, this new domain size is adequate as long as 
the overlap length, l(T m — T) > R(T,t w ), however, at 
temperatures just outside T < (T m — AT*) the typical 
domain size is determined by the overlap length 
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FIG. 8. x"(T) - x"ef(T) vs. T, measured on heating at 
0.25 K/min., after the sample has been cooled to 23 K at 
0.25 K/min. (circles) or 4 K/min. (squares) and kept at 23 
K for 3 600 s (solid symbols) or 30 000 s (open symbols). The 
pluses show a heating curve measured after cooling the sam- 
ple at 4 K/min. without a stop at 23 K. / = 1.7 Hz, h ac = 0.1 
Oe, Ag(ll at% Mn). 

l{T m - T) < R(T,t w ). (The overlap length, Z(AT), is 
a very rapidly decreasing function with AT). At lower 
temperatures, the cooling time again becomes the gov- 
erning parameter. However, just outside the overlap re- 
gion, T m — AT*, the actual cooling time, t ac , i.e. the time 
the sample has been kept within a region where the grow- 
ing domain size overlaps with the domain size achieved at 
temperatures just above, is of course considerably shorter 
than the cooling time, t c , of the continuous reference 
cooling procedure and thus R(T,t ac ) < R(T,t c ). During 
the continued cooling, t ac increases, but a temperature 
decrease at least of order AT* is needed to regain t ac « t c 
and a typical domain size equivalent to that of the con- 
tinuous cooling procedure. These arguments suggest, in 
accord with the experiments, that the disturbed x"(T)- 
curve at temperatures below T m first should increase to- 
wards and even cross the reference level as the temper- 
ature is decreased beyond AT*, go through a maximum 
deviation and then slowly approach x'ref further away 
from T m . Similar arguments as discussed above can also 
be applied to the heating curve at temperatures above 
T m , i.e. the dip should be succeeded by a temperature 
region outside T m + AT* where the relative susceptibility 
should become positive. The longer time the sample is 
kept at T m the deeper the dip and the larger the positive 
levels outside the overlap region will become. This con- 
clusion is further elucidated in an experiment where the 
sample is cooled, either with the rate 0.25 K/min. or 4 
K/min., to the temperature 23 K where it is subjected 
to different wait times before it is heated with the rate 
0.25 K/min. In Fig. 8 the ac-susceptibilities, measured 
on heating, are indicated by filled symbols after the wait 
time, t w = 3 600 s. Open symbols indicate the corre- 
sponding curves after a wait time, t w = 30 000 s. To 
separate the two cooling rates, squares are used to 
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FIG. 9. Relaxation rate S vs. logt, measured at 25 K af- 
ter having cooled the sample at 0.25 K/min. to 23 K and 
keeping it there a wait time tw before re-heating it to 25 K 
and applying the de-field after 10 s at constant temperature. 
t w = (filled squares), 3 600 s (open circles) and 30 000 s 
(solid circles). A reference curve measured after cooling the 
sample at 4 K/min. immediately to 25 K, and applying the 
field after 10 s is indicated by open squares. Hdc = 1 Oe, 
Ag(ll at% Mn). 

symbolize the faster one while circles are used for the 
slower. Also included in Fig. 8 is a reference curve ob- 
tained during heating after cooling to 20 K with the fast 
rate, 4 K/min. As can be seen, the system appears more 
reinitialized, at high enough temperatures, after a wait- 
ing time of 30 000 s. The effect is most pronounced in 
the case where a slow cooling rate was used. A wait time 
of 3 600 s is apparently too short realize a measurable 
reinitialisation. 



C. Relaxation rate experiments 

The usefulness of this fractal domain picture can be 
further examined by means of relaxation measurements 
where the observation time spans several decades. A 
zero-field-cooled relaxation experiment can be designed 
to give complementary information on the reinitialisa- 
tion observed in the ac-susceptibility curves of Fig. 7. In 
these measurements, as in the ac-experiments, the sam- 
ple is cooled with a rate of 0.25 K/min. to 23 K where 
it is subjected to a wait time, t w . The sample is then 
heated at the rate 0.25 K/min. but only to 25 K where 
the sample is kept at constant temperature and after a 
time of 10 s a magnetic field of 1 Oe is applied. The sub- 
sequent relaxation rate is shown in Fig. 9 for the three 
different wait times, t w = (filled squares), 3 600 s (open 
circles), and 30 000 s (filled circles). For reference pur- 
pose, the corresponding rate obtained after cooling the 
sample at 4 K/min. directly to 25 K and waiting there 10 
s before the field is applied, is included and indicated by 
open squares. It should also be mentioned that a similar 
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FIG. 10. Relaxation rate S vs. logi, measured at 31 K. 
The sample was cooled to 31 K, kept there a wait time 1 
000 s, then cooled to 29.7 K, kept there 3 000 s and finally 
re-heated to 31 K, where the dc-field was applied. Hdc = 1 
Oe, Ag(llat% Mn). 

measurement where the sample is directly cooled to 25 
K using a cooling rate of 0.25 K/min. gives a relaxation 
rate curve that is closely equal to the curve where a zero 
second wait time at 23 K is used (filled squares). The 
prime aging characteristic of the magnetic relaxation in 
spin-glasses is an inflection point in the relaxation rate 
at an observation time closely equal to the wait time 
before the magnetic field is applied, t w . A rather slow 
cooling rate, like 0.25 K/min., yields a broad and ill de- 
fined maximum in the relaxation rate at an observation 
time that may be called effective wait time or age. The 
three curves taken after a cooling rate of 0.25 K/min. 
in Fig. 9 show a continuous development towards the 
rapidly cooled (4K/min.) reference curve with increas- 
ing wait time at 23 K. In terms of the droplet model 
and in agreement with the discussions above, the behav- 
ior may be interpreted as follows; a longer wait time at 
T = 23 K should create equilibrated regions on larger and 
larger length scales, length scales that should approach 
the typical domain size achieved during cooling at higher 
temperatures. At 25 K, which on our experimental time 
scales is out of the overlap region, these regions looks ran- 
dom, i.e. with increasing wait time at 23 K, the system 
should appear more and more reinitialized at T = 25 K, 
just as Fig. 9 shows. 

The established way to show the possibility for a spin- 
glass to simultaneously carry several characteristic length 
scale, at one-temperature, is by temperature cycling 
experiments.!!!] The relaxation rate result presented in 
Fig. 10 is obtained by first cooling the sample to T\ = 31 
K where it is subjected to a wait time of 1 000 s. The 
temperature is thereafter changed to T2 = 29.7 K where 
it is kept for t W 2 — 3 000 s. When returning to T\ 
again after this temperature cycling, diluted regions of 
size R(T2, t W 2) have been reinitialized due to the domain 
wall movements at T2. The effect of a third wait time 
tu>3 = 10 s at Ti is a system which initially equilibrates 
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FIG. 11. a) Relaxation rate S vs. logt, measured at 27 K 
after a wait time of 3 000 s. Crosses, the sample was imme- 
diately cooled to 27 K, open circles, the sample was intermit- 
tently kept at 31 K for 3 000 s during cooling, b) Relaxation 
rate S vs. logt, measured at 31 K after having heated the 
sample back from 27 K and immediately applying the dc-field 
(open squares). Open circles shows a corresponding curve 
where the sample was kept at 23 K for 3 000 s instead of at 
27 K. Solid circles shows the relaxation rate after cooling the 
sample directly to 31 K and wait 3 000 s before applying the 
field. H dc = 1 Oe, Ag(llat% Mn). 

in these reinitialized regions up to length scales 
R(Ti,t w s). By choosing the temperature step and the 
wait times carefully it is possible to fulfill the condition 
R(Ti,t w3 ) < R(T 2 ,t W 2) < R(Ti,t w i). The two maxima 
in the relaxation rate in Fig. 10 then correspond to the 
non-equilibrium response due to the remaining two rel- 
evant characteristic length scales at T%, R(Ti,t w g) and 
R(Ti,t w i). 

It is also possible to show that several length scales 
can coexist in the spin-glass at different temperatures 
with relaxation measurements. Indicated by crosses in 
Fig. 11a is the relaxation rate for the temperature 27 K 
after a wait time of 3 000 s. Also shown in the same fig- 
ure, indicated by circles, is the relaxation rate after the 
sample has spent 3 000 s at 31 K before the aging at 27 
K. No difference in the measured relaxation rate can be 
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resolved between these two procedures. A length scale, 
caused by the wait time at 31 K, is however hidden in the 
second case. In order to reveal this, the sample is heated 
back to 31 K and immediately probed in a magnetic field. 
The subsequent relaxation rate is indicated by squares in 
Fig lib. Also shown by open circles in this figure, is the 
corresponding relaxation rate but with the second aging 
at 23 K instead of 27 K. As a reference, the relaxation 
rate probed directly after the initial wait time at 31 K 
is shown as filled circles. Basically, all three curves show 
an age of the system equal to 3 000 s only differing in 
the magnitude of the rate. This experiment, similar to 
the ac double memory experiment, shows that a spin- 
glass always recalls what has happened during cooling 
when it is re-heated to (or through) the temperature of 
an interrupted cooling process. 

V. CONCLUSIONS 

The non-equilibrium behavior of the dynamic suscep- 
tibility of a 3d spin-glass has been measured after specific 
thermal procedures which have been designed to expose 
consequences of a corresponding domain growth as pre- 
dicted in the Fisher and Huse droplet model of spin-glass 
dynamics.!!! A good qualitative agreement between these 
experimental results and the model predictions is found, 
which give support for the relevance of this real space 
model. However, there remain important questions to 
be answered. One crucial point within the model is to 
obtain an understanding of the difference between the 
initial fractal domain configuration after a quench and a 
corresponding snap shot of the equilibrated system with 
droplet excitations on all length scales. Is it a differ- 
ence of the fractal dimension of the instantaneous do- 
main walls in the two cases that signifies the difference? 
A prime remaining issue is also to theoretically establish 
if the concept of fractal domains on many length scales 



can be valid for a model Ising as well as for a real 3d 
spin-glass. 
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